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Based on the perturbation theory up to the second order, analytical asymptotic expressions for 
the variation of the population of hyperfine component /r = 0 particles in the ground state of spin-1 
condensates against a magnetic field B has been derived. The ranges of B in which the asymptotic 
expressions are applicable have been clarified via a comparison of the numerical results from the 
analytical expressions and from a diagonalization of the Hamiltonian in a complete spin-space. It 
was found that, For Rb, the two analytical expressions, one for a weak and the other one for a strong 
field, together cover the whole range of B from 0 to infinite. For Na, the analytical expressions are 
valid only if B is very weak or sufficiently strong. 
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1, Introduction 

The spinor condensates, as tunable systems with active 
spin-degrees of freedom, are rich in physics and promising 
in application. Since the pioneering experiment on spin-1 
condensates [l| the study of these systems becomes a hot 
topic. ii In general, in the study of spin-/ condensates, 
an important observable is the population of the particles 
lying in a given hyperfine-component /x = — / to /. These 
quantities are popularly measured in various experiments 
and are a key to relate experimental results to theories. [H 
Hence, the theoretical study of these populations is 
desirable. 

Let the hyperfine density of /x = 0 component in the 
ground state (g.s.) of spin-1 condensates be denoted as 
Po- When a magnetic field B is applied, the variation of 
po as a function of B is studied in this paper. A many- 
body theory, instead of the mean-field theory, is used so 
that the spin-degrees of freedom are treated rigorously. 
The emphasis is placed on the asymptotic behaviors of 
Po when H —>• 0 and —>■ oo. Based on the second order 
perturbation theory analytical asymptotic expressions for 
Po are derived. Numerical results from a diagonalization 
of the Hamiltonian are also obtained to compare with the 
results from the asymptotic expressions. In this way the 
effective ranges of B wherein the asymptotic expressions 
are applicable are clarified. Thus, within the effective 
ranges, po can be quantitatively and accurately known, 
and in general the variation of po as a function of B can 
be qualitatively understood. 


2, Hamiltonian in the spin-space 

We consider N spin-1 cold atoms trapped by an 
isotropic potential and subjected to a magnetic field. For 


spin-1 atoms the dipole-dipole interaction is extremely 
weak, and thus its effect on the coupling between the 
spatial- and spin-modes can be neglected (the quasi¬ 
spin-orbit coupling is not considered here). Let the to¬ 
tal interaction V = VqV 2 , Vq = “ ^j) 

and V 2 = C 2 J2i<j ~ ^j) where is the spin- 

operator of the i-th particle. [si Isf It was found that 
|c 2 /co| = 0.0046 for ®^Rb and 0.031 for ^®Na. Thus 
the spin-dependent force is nearly two order weaker than 
the central (spin-independent) force. Therefore, the spin- 
modes are much easier to get excited. Accordingly, the 
lowest lying levels would avoid spatial excitations but be 
dominated by spin-modes. With this in mind, it is as¬ 
sumed that the g.s. does not contain spatially excited 
modes, and all the particles fall into the same spatial 
state which is most advantageous to binding. This is the 
basic assumption of this paper. 

The Hamiltonian H is 


H = Ho + V + qB^fl 


( 1 ) 


where Hq — + ^^(g)) includes the kinetic 

and trap energies. Since the Hamiltonian conserves the 
total magnetization M, the linear Zeeman term has been 
neglected. The third term at the right of eq.(l) is for the 
quadratic Zeeman energy, where q is proportional to B^. 

When 5 = 0, not only M, the total spin S is also 
conserved. It has been proved that, for a spin-1 N — body 
system, the all-symmetric spin-state with a pair of good 
quantum numbers {S, M) is unique, where N — S must 
be even. Let the total spin-state be denoted as Usm- 
The set {'^sm} form a complete set for all-symmetric 
spin-states. Thus, under the above basic assumption, 
the g.s. can be written in the form as 111 


$gs — Tli(j)(vi)'dSgM 


( 2 ) 
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where 4>(vi) is the common spatial state for all the par¬ 
ticles, Sg = N for Rb, and Sg = M or {M -\- 1) for Na 
a N — M is even or (odd). In what follows M > 0 and 
N — M being even are assumed for convenience. 

When B ^ Q, S is not conserved but M is. Hence, 
different 'dsM having the same M but distinct in S are 
mixed up, and the g.s. can be in general written as 

$gs = ^ Cs'&SM = ni(()(ri)0gs (3) 

s 

Since the spatial degrees of freedom are considered as 
being frozen, the Schrodinger equation can be projected 
into the spin-space as 

J = 0 (4) 


Qn-2,N-2 

q[M^{2N+ 3)+ 2N^-7N^+ 5N-2] 

(21V- l)(2iV- 5) 

and Qn— 4 :,n = 0. 

On the other hand, for any total spin-state, the op- 

A A 

erator is equivalent to N — Nq, where Nq is the 

operator for the number of /r = 0 particles. Thus we 
have 


(0gs|H(|Qgs) 

q{Qgs\Qgs) 


{egs\N0\Qgs) 


(0gs|0gs) 


= N- Npo 


( 10 ) 


Inserting eq.® into (HOD, and making use of eq. dzD to 
Q, we have the hyperfine density with ^ = 0 as 


where the integration covers all the spatial degrees of 

a 2 

freedom. Making use of the fact that J2i<j ~ 

A 

N, where S is the operator of the total spin, the above 
equation becomes 

(iJ' - E')eg, = 0 (5) 

a 2 

Where TV' = TV' + iV' = ^Xc 2 S , H[ = q^Jg, 
E' = E-Nh- -f NXc 2 , X = f dr|^|4, and 

h = {(j)\ — -I- U{ri)\(j)). In what follows we work 

simply in the spin-space with the Hamiltonian El'. 


3, Asymptotic behavior of the hyperfine density when T? —>■ 0 


When T? —^ 0, TV( can be treated as a perturbation. 
Since Dsm is an eigenstate of TVq and the set {Dsm} 
is complete, the perturbative states can be expanded in 
terms of {iTsm}- A crucial point is the calculation of the 
matrix element (^Ts'M|TV(|■^TsM) = Qs'S- Making use of 
the fractional percentage coefficients by which the spin- 
state of a single particle can be extracted from iTsm,[ 12| 
the general expression of Qs'S can be obtained as given 
in eq.(2) to (7) of ref.fl^. It turns out that Qs'S can be 
nonzero only if S" = S' or S' ± 2. 

For Rb with C 2 < 0, the unperturbed g.s. in the spin- 
space is ^NM- Due to the limited choices of S' and S', 
the perturbative state up to the second order is 

0gs ~ ^NM + P''dN-2,M + P y t (6) 

A C2 (4iV — Dj 

where /?' = P = x^c^{ 2 N-i) ■ From 

the general expression of Qs'S^ 


Qn,n = qiM^ + TV" - N)/{2N - 1), 


(7) 


Qn-2,N = 


2N-1 


2((TV- 1)2 - - M^) 

(2TV - 3) 


(8) 


po ~ (TV2-M2)/[lV(2iV_i)] (11) 

— -^(f^P'QN-2,N + {P')'^{Qn-2,N-2 — Qnn)) 

where the terms higher than q^ are excluded. The first 
term is the g —^ 0 limit of po for Rb and is denoted as 
(/Oo)g=o, which has been found before as given in eq.(18) 
of the ref. 0. Thus eg. (1121) is a generalization of the 
previous finding for magnetic field zero to nonzero. It 
turns out that (TV^ — M^)/[TV(2TV — 1)] = (C^g^_^ 
where the Clebsch-Gordan coefficient is introduced. This 
is because, when q = 0, the total spin of the g.s. Sg = 
N. Thus, each particle must couple to a (TV — 1)—body 
system with total spin TV— 1. Thus the probability of each 
single particle in /j, = 0 is given by the above square. 

It is notable that (po)g=o is irrelevant to dynamics but 
depends decisively on M. In particular, (po)g=o is close 
to 1/2 when TIV is small, and is close to zero when M Ki N. 
Note that, for the frame upon it TIV is defined, M k, 0 
implies that Sg is lying close to the X-Y plane and there¬ 
fore nearly half of the particles are in /i = 0. While the 
case TVV ft! TV implies that Sg is lying close to the Z-axis 
and therefore the number of /r = 0 particles should be 
small. Note that the second term at the right of eq. m 
is positive (because /3 is negative when C 2 < 0). Thus, 
starting from (po)g=o, Po keeps increasing with B. This is 
natural because the g.s. would do its best to increase the 
number of p, — 0 particles so as to reduce the quadratic 
Zeeman energy. 

For Na with C 2 > 0, the unperturbed g.s. is Dmm- The 
related matrix elements of TV( are 

QMM = q{2MN + 2N + M)/{2M + 3) (12) 

^ 2q AM+ 1){N - M){N + M+ 3),i 

Q««,m = -1- 

(13) 

Qm+2,m+2 = (2M-b3)(2TlV-h7) 

[TH2(4TV + 2 ) + TVT(10TV - 5) -f lOTV - 6] 
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The perturbative state up to the second order is 


©gs ~ I^MM + l'^M+2,M — 7 


Qm+4,M+2 

Xc2(4M + 10) 


M+4,M 


where V — 'yi\ + Qm+2,m+2-Qmm n __ Qm+2,m 

wneie 7 —71^1-1- xc2{2N-i) I’ i — Xc2(2M+3)- 

Inserting ea. ([T5l) into ea. ifTOll . we have 


(15) 


Po 


(N - M)/[Ni2M + 3)] - J-{2j'Qm+2,m 

qN 

+ {'1')'^{Qm+2,M+2 — Qmm)) (16) 


where the terms higher than are excluded.. The first 
term is (p())„= n fo r Na. This form was first found in 2000 
in the ref. [l4 Il5| and also in [l^ derived in a different 
way later. Thus eq. (HU) is a generalization of the previous 
finding from B zero to nonzero. When M = 0 the g.s. 
is in a pure polar phase and every particle is in a singlet 
pair. In this case the first term becomes 1/3 as stated 
in ref. ip and in the Theorem I of ref. [l 6 | . Since the 
square of the Clebsch-Gordan coefficient for the singlet 
pair = 1/3) the probability of the particle in 

/j, is therefore 1/3. This is the physical background of this 
value. When M increases even a little from zero, the first 
term at the right of eq. (1161) will decrease remarkably from 
1/3 (say, when M increases from 0 to 2, (po)g=o decreases 
from 1/3 to ~ 1/7)- This remarkable decrease arises 
from the factor 2M + 3 in the denominator. It implies 
that the appearance of a few polarized particles among 
the big group of singlet pairs could cause serious effect. 
Thus, for the case with C 2 > 0, the sensitivity against M 
when M is small is notable. Starting from (po)g= 0 ) Po 
keeps increasing with B (because 7 is negative) as in the 
previous case. 


4, Asymptotic behavior of the hyperfine density when B —^ 00 


When B —>■ 00 , Hq, rather than H[, can be treated 
as a perturbation. Since the Fock-state |fVi, Nq, fV_i), in 
which Nfi particles are in the /r-component, is an eigen¬ 
state of H[, the perturbative states can be expanded in 
terms of them. Since N±i = {N — Nq± M)/2, the Fock- 
state can be simply denoted as |iVo) when M is fixed. Let 
(A^Q|iLo|iVo) = P/v^Wo- Making use of the formula, say, 
given in eq.(A5) of ref. 17|, we have 


PnoNo = ^{M^ + N + No + 2NNo - 2Nl) (17) 


Pno-2,No = ^[fVo(lVo-I)((fV-fVo+2)'-M2)]5 (ig) 


Pno+2,No = ^[{No + l)(iVo + 2){{N - Nof - M^)]^ 

(19) 


Otherwise, P^'Afo = 0. 

Note that, when C 2 —>■ 0, the number of /r = 0 parti¬ 
cles in the g.s. will be maximized so as to minimize the 
Zeeman energy. Thus, under the conservation of M, the 
leading term of Qgs should be |iV —M). Accordingly, the 
perturbative state up to the second order is 


Qgs^ \N-M) + 5'\N-M-2) 

+ ^ ^A^-M-4,X-M-2 

where < 5 ' = < 5(1 + p^-m-^.^-m- 2 -Pn-m.n-m 

Pn-m-2,n-m 

2q 

Inserting Qgs into eq. (HI, we have 


( 20 ) 
<5 = 


Po 


N-M 
N 



( 21 ) 


where higher order terms have been neglected. This for¬ 
mula holds for both Rb and Na. Since the second term 
at the right of eq. (HU is negative, the first term appears 
as the upper bound. 


5, Applicability of the asymptotic expressions for C 2 


In order to clarify the applicability of the asymptotic 
expression, we have to find out the exact solutions of H'. 
It is reminded that the set of eigenstates of Hq, {'dsM}) 
is complete for all the symmetric spin-states, therefore 
they can serve as the basis functions for the diagonaliza- 
tion of H' in the spin-space. The related matrix element 
{'&S'm\P['\'&sm) = + I) + Qs'S) where Qs's 

is given in the ref.[l3|, and X = f c?r|(/|^ depends on 
the interaction and the trap. To obtain numerical re¬ 
sults, ?7(r) = is assumed. Note that X does 

not appear in the leading terms of all the asymptotic ex¬ 
pressions, and we consider only the cases with a large N. 
Thus, it is reasonable to evaluate X under the Thomas- 
Fermi approximation (TFA). For details, we refer the 
reader to ref. 18|. After the diagonalization one can ex¬ 
tract Po from the eigenstates of H' as given in eq.(9) 
of ref. [ 1 ^. Numerical results of po from the asymptotic 
expressions and from the diagonalization are compared 
below. 

It turns out that, when q = 0 and q = 00 , the resultant 
Po from the perturbative approach and from the exact 
diagonalization of H' are identical. Furthermore, we al¬ 
ways have (po)g->oo > po > (po)g=o (where the equality 
holds only if M = TV). Thus the variation of po versus q is 
strictly restricted in a domain which can be known in ad¬ 
vance, and Po keeps increasing with q inside the domain. 
The increasing arises because the g.s. would like to have 
more p = 0 particles to reduce the quadratic Zeeman en¬ 
ergy. We introduce a ratio a so that M = aN. For Rb, 
we know from eq. dn and that (po)g=o « 5 ( 1 -^^) 
and (po)g^oo = 1 — 0 !. In what follows, a is given at 1/2 
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and 0.001, respectively. Accordingly, the magnetization 
is half and nearly zero. The case a ~ 1 (corresponding 
to a nearly full magnetization) is trivial because both the 
lower and upper bounds are close to zero, and nearly no 
fj, = 0 particles will emerge. 

Examples of the variation of the densities po versus 
B with N = 10000 and a = 1/2 (0.001) is shown in 
Fig.l (Fig.2). Accordingly, po is increasing from 3/8 to 
1/2 and from 1/2 to 0.999, respectively. To measure the 
deviation between the exact and asymptotic results, we 
define Xj = \{po)exact-iPo)asym,i\ where {po)exact denotes 
that the density is from the exact diagonalization of H/ 
while {po)asym,i Is from the asymptotic expression ea. lfT^ 
(if i = 1) or from ea. d^ (if i = 2). To give a quantitative 
description, we define Bi at which Xi = 0 . 01 , and B 2 at 
which X2 = 0.01. We found that, when B < Bi {> B 2 ), 
xi (X 2 ) is even smaller than 0.01. Thus, we can say that 
the effective range of in which eq. uni) is very close to 
be exact, is [ 0 , i?i], and the effective range for ea. d^ is 
[.62,00]. 

The case with a half-magnetization (a = 1/2) is shown 
in Fig.l. In Fig.la ( 61 , 52 ) = (SdSmG, 361mG). Thus, 
the effective range for (po)asym.i and the effective range 
for {po)asym,2 together cover nearly the whole range of 6 . 
It is emphasized that both the lower and upper bounds 
depend only on a but not on N and/or w. Therefore, 
when N and/or uj are changed, the variation of the curves 
is limited because they are fixed at their two ends and 
they must keep increasing with B. Nonetheless, 61 and 
B 2 will therefore be changed. For an example, the uj in 
la is two times the lo in lb. Accordingly, the curves in 
la as a whole shift to the right. Numerical examples of 
( 61 , 52 ) with different N and ut are given in Table I. 

We know from the Table that B 2 is close to 5i in gen¬ 
eral, and therefore the asymptotic expressions are nearly 
valid in the whole range of 5. They will both be larger 
when N and/or uj are larger, or vice versa. Nonetheless, 
they are more sensitive to w rather than N. When w is 
larger, the spin-texture of the g.s. would have a stronger 
ability to resist the field. This causes a shift of the curves 
to the right. 

The case with a nearly zero magnetization (a = 0.001) 
is shown in Fig.2. It was found in 2a (2b) that both 5i 
and B 2 are very close to a critical point Bcru = 446mG 
(292 toG) at which (po)exact undergoes a sharp change. 


TABLE I: Bi and B 2 in mG for Rb with a = M/N = 1/2. 
Two cases of N and three cases of ui are considered, where 
LUo = 300 X 2-k/s. 


The change is so sharp that the second order deriva¬ 
tive of (po)exact against 5 tends to 00 when N ^ 00. 
Obviously, this implies a phase transition. It turns out 
that the two curves {po)asym,^ and ipo)asym,i intercept 
when 5 ^ Bcrtt, and this is a common feature for the 
g.s. with C 2 < 0 and with a nearly zero-magnetization. 
Thus the effective ranges for the two asymptotic expres¬ 
sions together cover the whole range of 5, and the two 
expressions together provide a perfect description of po- 
Numerical examples of Bcrtt are given in Table II. 

Similar to the previous case, a larger N and/or a larger 
w lead to a larger Bcrit, and the vice versa. When 
5 > Bcrit, the curve of {pQ)exact becomes horizontal. 
It implies that, after the phase transition, the system ar¬ 
rives at its eventual status (this status is described by 
the Fock-state |M, A^ — M, 0), where N — M particles 
have p = 0 while the rest have p = 1). In this status 
the number of p = 0 particles has been maximized, and 
therefore no more change is allowed when 5 increases fur¬ 
ther. Note that the sharp change will become ambiguous 
when M is larger. Comparing 2b with 2a, it is clear that 
Bcrit will shift to the right when uo increases. 


6, Applicability of the asymptotic expressions for C 2 > 0 

Let us predict two features of spin-1 condensates with 
C2 > 0. 

(i) When 5 = 0, the total spin-state of the g.s. has 

S = M and is denoted as i9mm oc Sx^ , 

where xi denotes the spin-state of a single particle with 

p = 1, {xx)o denotes a singlet pair, and 5 is a sym- 
metrizer. If there are another forms, they are identical 
due to the uniqueness of the eigenstate. [lOj Thus, the g.s. 
is a mixture of a group of polarized particles together 
with a group of singlet pairs. Obviously, the M particles 
in the first group are stable against 5, while the N — M 
particles in the second group are not. Note that every 
singlet pair is situated under the same environment, thus 
they have similar ability to resist 5. It is possible that 
all the pairs might begin to be broken when 5 increases 
and exceeds a certain value, and therefore a sharp change 
of the spin-texture will occur. 

(ii) In general, the stability of the g.s. depends on the 
gap (the energy difference between the g.s. and the first 


TABLE II: Bcrit in mG for Rb with a = 0.001. In this case 
5i ~ 62 ~ Bcrit- uJo is referred to Table 1. 



N=10000 

N=20000 


N= 

10000 N=20000 

uJo/2 

226, 238 

260, 274 


292 

337 

OJo 

343, 361 

393, 414 

OJo 

446 

514 


520, 548 

596, 629 

2(Uo 

678 

778 
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excited state). When B = 0 the g.s. has S = M while 
the first excited state has S' = M + 2, therefore the gap 
is oc (4M + 6) for C 2 > 0, whereas this factor would be 
(47V — 2) for C 2 < 0. Therefore, the gap for Na is much 
smaller than that for Rb when a is small. In this case, 
the g.s. of Na is highly unstable and the feature of the 
system will depend on a sensitively. In other words, the 
singlet pairs will become more fragile when a —>■ 0 . 

We found from eq. (HU) and ea. (l2T]) that the lower and 
upper bounds are (po)(j=o = (1 ~ a)/{2aN + 3) and 
(po)q-K>3 = 1 — a, respectively, and they are identical 
to the numerical values from the exact diagonalization 
of iJ'. In particular, when a = 0, (po)<?=o = 1/3 and 
(po)ij=oo = 1- Note that, when a = 0, every particle 
is in a singlet pair. Recalled that the Clebsch-Gordan 
coefficient Therefore, when 

a particle is in a singlet pair, the probability in p, is 
1/3. This leads to (po)q=o = 1/3. On the other hand, 
(po)q=oo = 1 implies obviously that all particles are in the 
p = 0 component. This results from the minimization of 
the quadratic Zeeman energy. 

For any a, we have (po)q=o < Po < (po)q^oo, and 
po keeps increasing in between as before. Note that the 
(po)q-i-oo of the two species Rb and Na are the same, but 
the (po)q=o of Na is always lower than that of Rb. Thus, 
when q varies, the po of Na varies in a broader domain. 
Note that the a—sensitivity of (po)(j=o is embodied in 
the factor 2aN which appears in the denominator (Say, 
when N = 10000, and a = 0 and 0.001, respectively, 
(po)(j=o = 1/3 and 1/23. Thus a small change in a leads 
to a big change in (^ 0 ) 5 = 0 )- 

For the case of half-magnetization (a = 1/2), the vari¬ 
ation of Po versus B is shown in Fig.3 where N = 10000 
and u) = ujo = 300 x 27r/s. There is a sudden uprising 
in (po)exact as mentioned in point (i) taking place when 
B = Bcrit = 114mG. It turns out that {po)asym,i is iden¬ 
tical or extremely close to (po)exact when B < Bcrit, but 
deviates from {po)exact rapidly when B > Bcrit- There¬ 
fore, in this case, we have Bi Ri Bcrit and therefore the 
effective range for (po)asym,i is {0, Bcrit)- Numerical ex¬ 
amples of Bcrit and B 2 (the latter is much larger) are 
shown in Table III. 

Thus, both Bi and B 2 will be larger when N and/or w 
are larger as before, or the vice versa. In fact, Bcrit mea¬ 
sures the ability of the singlet pairs to keep themselves 
unbroken against B. It is obvious that this ability will 


TABLE III: Bi ~ Bcrit and B 2 in mG for Na with a = 1/2. 
Refer to Table I. 


become stronger when the trap is stronger (w is larger). 
It is interesting to see that this ability will also become 
a little stronger when N increases. 

Recalled that the Bcrit for Rb marks the maximization 
of the p = 0 particles. Now, the Bcrit for Na marks the 
solidity of the singlet pairs. This explains why the curve 
is horizontal when B > Bcrit for Rb, and when B < Bcrit 
for Na. In the latter case the pairs are kept unbroken and 
the spin-texture remains unchanged. 

The case of nearly zero-magnetization {a = 0.001) is 
shown in Fig.4. It is shown in Fig.4a (for B < ImG) 
that {po) exact increases very fast when B is ranged from 
0.2mG to ImG. One can define a critical strength Bcrit 
at which the second order derivative of {po)exact against 
B arrives at its maximum. In Fig.4a Bcrit = 0.26toG. 
It is further found that the effective range of B for 
{po)asym,i to be applicable is from 0 to Ri = 0.29 toG. 
Thus Bi is again very close to Bcrit- Once B > Bcrit, 
although both {po)asym,i and {po)exact increase sharply, 
they deviate more and more from each other. For the 
case that a is very small, the effective range for {po)asym,i 
is very narrow. As an example, if a is further reduced 
from 0.001 to zero, Bi would be reduced from 0.29mG to 
0.04mG, and Bcrit is reduced from 0.26mG to zero. Thus, 
the singlet pairs will become very fragile when a —^ 0 as 
predicted in (ii). 

How Bi and B 2 are affected by N and uj is shown in 
Table IV. 

Table IV, Bi and B 2 in mG for Na with a = 0.001. 
Refer to Table I. 


7, Summary 

When B tends to zero and infinite, the asymptotic 
forms of the hyperfine density po against B denoted as 
{po)asym,i and {po)asym,2, respectively, has been obtained 
analytically based on the second order perturbation the¬ 
ory. Numerical calculation via an exact diagonalization 
of H' for the density denoted as (po)exact has been per¬ 
formed to clarify the effective ranges of B wherein the 
asymptotic forms {po)asym,i are applicable. The two 
ranges for i = 1 and i = 2, respectively, are specified 
as (0,i?i) and {B 2 ,oo). It turns out that the two limits 
(po)g=o and (po)g=oo given by the analytical expressions 
are identical to those from the exact diagonalization of 


TABLE IV: Bi and B 2 in mG for Na with a = 0.001. Refer 
to Table 1. 


N=10000 N=20000 
uJc/2 77, 226 90, 259 

Wo 118, 342 136, 393 
2wo 176, 519 204, 595 


10000 20000 
uJo/2 0.17, 163 0.18, 187 

Wo 0.26, 248 0.27, 282 

2wo 0.39, 375 0.41, 428 
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H'. They together provide the lower and upper bounds 
for po) and po is monotonously increasing with q between 
them. Thus, in any case, we can have a rough impression 
on po- 

For Rb, both Bi and B 2 are in the order of lO^mG, 
and they are close to each other when the parameters 
vary in a broad domain that are accessed frequently in 
experiments. Thus, po can be accurately known when 
B lies inside one of the ranges, or can be roughly known 
when B lies between Bi and B^. In particular, when a = 
M/N is small, (i) A critical strength Bern is found which 
marks the realization of the eventual status, in which the 
number 0 / p = 0 particles is maximized, (ii) Bi, B 2 , 
and Bcrit are close to each other, and therefore the two 
asymptotic forms together cover nearly the whole range. 

For Na, a critical strength Bcrit is also found but has 
another implication, it marks the sudden breaking of the 
singlet pairs. At a weak held po is highly sensitive to 
a when a is small. Note that, when a increases from 
zero, a few polarized particles will emerge among the nu¬ 
merous singlet pairs and lead to a remarkable increase 
of Bcrit- It implies that the singlet pairs will become 
more solid thereby. Thus the emergence of a few specihe 
particles can cause serious effect, similar to the serious ef¬ 
fect of a few impurity appearing in well organized crystal 
structure. The underlying physics deserves to be further 
studied. 

In general, when the spin-texture undergoes a sharp 
change (this appears as a phase transition when N —>■ 
00 ), {po) exact is able to describe this change. However, 
the perturbation theory fails to describe such a sharp 
change. This is the main shortcoming of the perturba¬ 
tion theory and is the reason that {po)asym,i will begin 
to deviate sharply from {po)exact when q is close to the 
critical points. 
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FIG. 1: (Color online) The three densities {po)exact (solid), 
(po)asi;m.i(dash), and (po)as!/m. 2 (dash-dot) of Rb versus 
BjmG. M = N/2. The trap is assumed to be 
ipo)asym,i is effective (close to {po)exact) when B lies in 
(0,Bi), while {po)aBym ,2 is effective in {B 2 ,oo). Bi and B 2 
are marked by two vertical dotted lines. Note that the upper 
panel has a larger u>, and accordingly the curves shift to the 
right. It implies that, for a stronger trap, a stronger B is 
needed to raise up the number of /r = 0 particles. 
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FIG. 2: (Color online) The same as Fig.l but with M = 
A'^/1000. There is a sharp change in {po)exact at Bcrit marked 
by a vertical dotted line. B\ « Bcrit ~ B 2 is found. 



FIG. 3: (Color online) The three densities of Na versus 
BfmG. Refer to Fig.l. 
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. (b) Na / 

N=10'^, VI = N/1000 / 
co = 300x27i/s 
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FIG. 4: (Color online) The same as Fig.3 but with M = 
A^/1000. The range of B is 0 —^ ImG in 3a where is 
marked, whereas B is 1 —^ 300mG in 4b where B 2 is marked. 
ipo)asym ,2 doBS not appear in 4a, while {po)aaym,i does not 
appear in 4b. 








